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In metallic ferromagnets, an electric current is accompanied by a flux of 
angular momentum, also called spin current. In multilayers, spatial variations 
of the spin current correspond to drive torques exerted on a magnetic layer. 
These torques result in spin precession above a certain current threshold. The 
usual kind of spin current is associated with translation of the spin-up and spin- 
down Fermi surfaces in momentum space. We discuss a different kind of spin 
current, associated with expansion and contraction of the Fermi surfaces. It is 
more nonlocal in nature, and may exist even in locations where the electrical 
current density is zero. It is larger than the usual spin current, in a ratio of 10 or 
100, at least in the case of one-dimensional current flow. The new spin current 
is proportional to the difference A/I ~ 10~ 3 eV between spin-up and spin-down 
Fermi levels, averaged over the entire Fermi surface. Conduction processes, spin 
relaxation, and spin-wave emission in the multilayer can be described by an 
equivalent electrical circuit resembling an unbalanced dc Wheatstone bridge. 
And A/I corresponds to the output voltage of the bridge. 
I. INTRODUCTION 

In magnetic multilayers, the s-d exchange interaction causes 1 an increase of 
the Gilbert spin-damping parameter a near an interface. In the presence of a dc 
current normal to layers, it also generates 1-5 drive torques on magnetic spins 
near the interface; these lead to a spin precession above a certain current thresh- 
old. This current-induced precession has been observed recently 6 ~ 8 in Co/Cu 
multilayers. Such device, which we call a SWASER, is the magnetic analog of the 
semiconductor laser diode. Current-induced switching between opposite static 
magnetic-spin directions has also been observed 7 ' 8 , at low magnetic fields; it 
happens probably through a temporary current-induced precession. Microwave- 
frequency ac currents have also been used 9 to excite precession in a multilayer. 
Detailed numerical as well as analytical calculations of dc current-induced pre- 
cession and switching have been done 10 for various thin-film geometries and 
anisotropy-field cases. A physically different mechanism has been proposed 11 
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for current-induced switching. 

The layer configuration most common in theory and experiments is shown 
in Fig. la. Magnetic layer Fi prepares the spin of conduction electrons in a 
fixed specific direction. It must be 1 thicker than a spin-diffusion length, in 
order to ensure a ratio a\ of spin-up to spin-down current as different from 
one as possible. The non-magnetic spacer N must be thinner than a local 
spin-diffusion length, to preserve spin memory for electrons crossing it. The 
precessing magnetic layer F 2 must be thin, as the current-induced drive torque 
is only applied 2 to a region of F2 of thickness ir/(k^ — fcj.) near the N — F2 
interface. Here, fcj and k[ are the spin-up and spin-down Fermi wavenumbers 
in F2, both assumed real- valued. Thick non-magnetic layer N2 acts as a current- 
return lead. We introduce 1 spatial coordinates x,y,z, with x normal to layers, 
and x=0 at the N — F 2 interface (Fig. la). We also use spin coordinates 

^2) 2/2, Z2, with Z2 parallel to the precessing atomic spins S2 in F2, and X2 in the 
plane of Z2 and of the precession axis (Fig. lb). The direction of the precession 
axis, which is also the direction of the total effective magnetic field, is assumed 
parallel to the fixed atomic spins S\ in F\. The Si direction is arbitrary. The 
precession-cone angle between Si and S2 is called 9. 

The current-induced torque f on the spins of layer F2 , responsible for main- 
taining the precession, can be calculated 4 from the angular-momentum current 
P, also called spin current. The important components are along X2- 

T X 2 = P X 2{X = 0) - P X 2(X = 00). (1) 

Here, P X 2(x = 0) is calculated in F 2 at the N — F 2 interface, and P X 2(x = 00) 
at a point of F 2 to the right of the region < x < 7r/(fcj — fcjj where most of 
the torque is applied. 

The traditional spin current arises from Fermi-surface translation, and is 
proportional to the local electrical-current density. The main purpose of the 
present paper is to describe a different kind of spin current, which comes from 
Fermi-surface expansion or contraction, and which can exist even where the 
current density is zero; this new spin current turns out to be much larger than 
the traditional one, at least when the current flow is one-dimensional. Another 
purpose is to demonstrate the equivalence of the Slonczewski and SWASER 
theories. 

II. SPIN CURRENT FOR ONE ELECTRON 

First, we consider free-electron states of wavevector fcjy in N which are eigen- 
states of the z 2 component of spin. If k% > 0, they enter F 2 from N, and their 
form in F2 is 



exp{ik^f) 2k 



(2) 
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The wavevectors fcj , ki are such that ip-\ , tpi have the same energy. These 
states are assumed to have most of their norm in the layer N of volume Vn- 
Momentum conservation along y and z at the iV — F 2 interface gives ky = k^ = 

and = k\ = k\. We neglect reflections at the F 2 — N2 interface. 

A part of the amplitude of the actual spin-up and spin-down waves tp- 
entering F2 from N with k x > originates 1 in F\. Hence, the spin density 
< s5(f — r ) > for that part is parallel or antiparallcl to S\, and at an angle 
9 to the Z2 axis. These states are linear combinations of - 0T)V'|- With z 2 as a 
spin quantization axis as before, we can write in F 2 : 



ft(r)cos(9/2) 
-/ i (r)sm(^/2) 



/;(f)cOs(0/2) 



(3) 



We assume 9^1 radian, so that ip + ,ip- can still be called spin- up and 
spin-down. The quantities T+ = Ak\k^ /(fcf + fc|) 2 ,T_ = 4fc^fcf '/(fcf + k [ x ) 2 
are the transmission coefficients at the Fi—N interface, assuming F\ to be made 
of the same material as F 2 . These quantities are present because we consider 
only the parts of the waves in F2 originating in F\. To simplify the theory, we 
approximate T + ~ T_ ~ T. 

For any state i[>, Slonczewski gives 4 an expression for the components 
Px2,P y 2 of the spin current: 



a(r) 
b(r) 



P 
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Using this on the k x > states ip + ,ip- of Eqs.(3), we obtain in F 2 , respec- 
tively: 



(p X 2(x ))+ = A sin envUvm^fcos{ki-ki) X 

4 V N (k? + ki)(k? + k l x ) 

(PA*))- = ~{Px2{x)) + . 

Here, v\ and v^. are electron- velocity components in F2. The cosine factor 
comes from the precession 1,A of the electron spin around the s-d exchange field 
set up by S2, as the electron moves toward increasing x in i 7 ^. We assume 
kjy < k^,ki so that the k x > states are not totally reflected at the N — F 2 
interface. 

On the other hand, the k^ < states originated in N 2 before traversing 
F 2 toward left. Because of the absence of a magnetic layer like Fx to polarize 
their spin along a specific fixed direction oblique to z 2 , we have P x2 = for 
these states, at any location in F 2 . The same applies to the part of the states 
in F2 that have k x > after reflection at the N — F2 interface. Also to the 
part (ignored earlier) that has k x > after reflection at the F\ — N interface, 
if T+~T_. 

In thermal equilibrium, k x > spin-up and spin-down states ip+, ip- of the 
same energy are equally populated, and Eqs. (5) show that their P x2 cancel 



3 



each other in pairs. The total P x i of the electron distribution vanishes, as 
expected from the second law since P X 2 is related to damping and drive torques. 
This cancellation does not happen in general for component P z2 , which plays 
no important role anyway. 

For the same reason, when out of equilibrium, only states near the Fermi 
level contribute to the total P X 2- 

III. TOTAL SPIN CURRENT IN CASE OF 
TRANSLATED FERMI SURFACES 

The traditional calculation 4 of the total spin current assumes a non-zero 
local current density j x , corresponding to translations of the spin- up and spin- 
down Fermi surfaces in N along the k% axis, by amounts <5fc+, 5k~ (Fig. 2a). 

Because of the cancellation mentioned above, only states located between 
the spin- up and spin-down Fermi surfaces in N contribute to the total P X 2, at 
zero temperature. Also, the blurring of the electron distribution near the Fermi 
level, happening at non-zero temperature, turns out to have no effect on the 
total P X 2- 

We integrate the P X 2 of Eqs. (5) over all number elements d 2 N of electron 
states, from a = to a = ir/2, corresponding to the k^ > half of the Fermi 
surface in N. For given energy, kl,k x are functions of a. As explained in the 
last section, only the states on that half have non-zero P X 2 ■ We call "EP X 2 the 
total spin current: 



„.„ , . Ti 2 sin6 , , , . m (2fci v ) 2 (fc| + k^)cos({k\ - k!i)x)cosasinada 

T,P x2 (x) = - — 9 kt(Sk+-dk~) / T y ^ ^ ^ . 

167r 2 m Ny X 'J {k N + k V ){k N + k i x) 

(6) 

For any given x 3> 7r/(fcf— fcjj, the argument of cos{{kl l . — k x \)x) in the integral 
varies rapidly with a. Hence, the average of this cosine over half a Fermi surface 
is nearly zero. In turn, this means SP^W — when x 3> 7r/(fcj — k{). Thus, 
the last term of Eq. (1) vanishes. Only Y.P x2 (x = 0) matters. 

An alternate, but equivalent, approach 1 to the problem of current-induced 
torques replaces the spin current by the spin imbalance, defined as the difference 
A[i = jU-f — m ~ 1CT 4 eV between the spin-up and spin-down Fermi levels at a 
given point of the Fermi surface. As before, only the electrons on the k^ > 
half matter, because some of them have their spin along S\ , and because they 
enter F 2 . 

When Afi < 0, the quantity A/i "pumps" the electron gas, just like the 
difference between valence-band and conduction-band Fermi levels does 12 in a 
laser diode, resulting in stimulated emission of spin waves, i.e., in precession. 
The current threshold where precession begins is given 1,6 by Afi+huj = 0, where 
lu is the angular frequency of precession, at any temperature. A mathematically 
identical condition applies 12 to the threshold of a laser diode. 

In the case of Fermi-surface translation, A/x is proportional to <5fc+ — 5k~ or 
to j x : 
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Here, ct\ = jl/j x is the local ratio of spin-up to spin-down current densities 
in N, and the electron density per unit volume in N. We average k x over 
the k x > half of the Fermi surface, obtaining k x ~ fcjv/2. Similarly, T av is 
some average of T over that half. As before, the T av /2 factor is present because 
only transmitted k^ > waves are active for SWASER action. We combine 
Eqs. (6,7) and obtain 



^Px2(x = 0) k^sinO 2 f n ^ 2 (2k x ) 2 (kl + k^.) cosa sina da . . 

J / t at . ttttt^t . ; in ■ (°) 



A/i " 8tt 2 T at ,7 (fc^ + fcl)(^ + fci) 

We see that £P X 2 and A/i vary in a constant ratio when j x or ai = jl/j x 
changes. This reflects the fact that the Slonczewski theory, based on Y,P x2 , 
is nearly equivalent to the SWASER theory based on A/z. The only significant 
difference is that the first theory uses only bulk damping, and the second theory 
only surface damping. 

IV. TOTAL SPIN CURRENT IN CASE OF EXPANDED 
OR CONTRACTED FERMI SURFACE 

There exists a different kind of spin imbalance, associated with isotropic 
expansion or contraction of the Fermi surface rather than translation (Fig 2b). 
It is the one originally introduced by Aronov 13 and by Johnson and Silsbee 14 . 
This quantity, which we denote 3 by ATI, is an average of fi-\ — U[ over the whole 
Fermi surface, not just one half. Since it depends on j\ and j\. values over a 
range of x of the order of a spin-diffusion length, we call it the nonlocal spin 
imbalance. It may exist even at locations where the electrical current density 
is zero. Like A/i < 0, Ajl < can be used 3,6 to induce spin precession in a 
multilayer. For the configuration of Fig. la, ATI has the same sign as Afi. 

We can calculate 3 ATI for the multilayer of Fig. la at the N — F 2 interface 
by solving the spin-diffusion equations in one dimension, assuming the F± and 
N2 thicknesses to be much larger than a local spin-diffusion length l sr , and the 
N, F 2 thicknesses much smaller than a local l sr . A formula for ATI was given as 
Eq. (5) of Rcf. 3. By comparing the formula to our Eq. (7) for A/i, one obtains 
roughly 

< (^±) F1 + (^±) N2 . (9) 
A/i l sr l sr 

Here, A^ and l sr are the spin-down mean free path and the spin-diffusion 
length, and Fl, N2 refer to layers Fi and iV 2 . We see from Eq. (9) that the 
ATI mechanism is usually dominant over the A/x mechanism, since we have 
KJl sr ~ 0.1 — 0.01 in most magnetic materials 15 . This is helped by the fact 
that Aj. < A T in Co, Ni, Ni-Fe. 
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The spin-current component P x2 of one electron is still given by Eq. (5). 
It is still true that only electrons located between the spin-up and spin-down 
Fermi surfaces (Fig. 2b) contribute to the total spin current of the electron dis- 
tribution, which we denote now by to differentiate it from the translation 
case. We consider a momentum-space volume element d 2 k at an angle a to the 
fcjy axis, containing d 2 N states of one spin: 

d 2 k = k% 2tt sina da Jfcjv (10) 
d 2 N = V N d 2 k/(8n 3 ). 

Here, 5k^ is the radial distance between Fermi surfaces, now independent 
of a. It is related to Ajl by 

Ajl=^-k N Sk N ^. (11) 

The quantity T av is a slightly different average of T than T av . Using Eqs. (10, 
11), we integrate the P x2 of Eqs. (5) over all number elements d 2 N to obtain 
YjP X 2- As before, only the k^ > half of the Fermi surface contributes: 

T,P x2 (x = 0) _ k N sin6 2 r /2 T (2k^) 2 (kl + k\.) sina da 
AMz = 0) ~~ 16tt 2 T^Jo (few + kl)(k« +kj) 

This formula for Fcrmi-surfacc expansion or contraction differs from that of 
Eq. (8) for translation by an extra 1/2 factor in the front, as well as by a missing 
cosa in the integral. Therefore, the ratios of spin current to spin imbalance given 
by the two equations are roughly the same. Just as Ajl dominates over A/z, 
^P X 2 is dominant over T,P X 2 if one-dimensional flow is assumed. This implies 
a considerable reduction of the threshold current where precession starts. 

As before, we have SP x2 ~ for x 7r/(fcf — fcjj. Only T,P x2 (x = 0) 
matters in Eq. (1). 

Most multilayer samples used in recent SWASER experiments 6-8 have cur- 
rent leads flaring out immediately to much wider effective cross section area 
than that of the multilayer itself. Thus, current flow is not one-dimensional. 
Because of the increased volume available for spin relaxation in these leads, A/J 
and T,P x2 are considerably reduced below their value in one dimension assumed 
in the present paper, for given current density in layer F 2 . Also, layer F\ is too 
thin, leading to the same result. On the other hand, A/i and YjP x2 in F 2 are 
not affected by these factors. 

In order to maximize \A~jl\ and IXP^I and to verify the predictions of the 
present paper, the cross-section area of layer N 2 , which plays the role of a lead 
(Fig. la), should remain equal to that of layer F 2 , over distances from layer F 2 
at least equal to a local spin-diffusion length l sr . And the F\ thickness should 
be at least a local l sr , again with constant cross section area. Like us, Fert and 
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Lee 16 find that A/i is not determined by interface scattering in the multilayer 
itself. 

The importance of Ap and T,P x2 can be compared to that of A/i and T,P x2 by 
measuring the threshold current of such one-dimensional sample, and comparing 
to that of a traditional sample with more multidimensional flow. Another way 
would be to realize the so-called zero-current SWASER 3 , which is based on ATI 
and T,P x2 alone. 

V. EQUIVALENT ELECTRICAL CIRCUIT 

Electrical conduction in the multilayer of Fig. la can be described by an 
equivalent dc electrical circuit (Fig. 3). It may be used to calculate Ajl. 

The upper horizontal resistors represent conduction in the spin-down band 
in F\ and N 2 , including both interface and bulk scattering, and the lower hori- 
zontal resistors conduction in the spin-up band. The vertical resistors represent 
electron-spin relaxation between the two bands. The dc electrical motor repre- 
sents the interaction between electrons and spin waves in F 2 . The dc current 
through the motor corresponds to the rate of electron spin flip in F 2 associated 
1 with spin- wave generation, i.e., with inducing spin precession. And A~p/e at 
the N — F 2 interface is represented by the dc voltage between points D and U. 
The counter-electromotive force of the motor is equal to Uu/e, and its internal 
resistance, shown as Rm, is inversely proportional to the number n m of spin- 
wave quanta (magnons) in F 2 . This circuit (Fig. 3) can be justified on the basis 
of the spin-diffusion equations 3 and of the formula 1 for the rate of production 
of magnons in F 2 . 

Consideration of the equivalent circuit, like that of the spin-diffusion equa- 
tions 3 , shows that the SWASER is affected by conduction processes in Fi and 
N 2 over a range of x of at least one spin-diffusion length, in one dimension. It 
also shows that A// and T,P x2 could be maximized, for given j x , in the case of 
Fermi-surface translation by having a very short spin-relaxation time in layer 
N 2 as compared to other layers. On the other hand, a long relaxation time in all 
layers is best in the case of expansion or contraction. This (large) dependence 
of YiP x2 on spin relaxation in the various layers is absent from Eqs. (31-33) 
of Waintal et al. 17 , since they ignore spin relaxation and the spin-diffusion 
equations. Also, they ignore N 2 . 

VI. CONCLUSIONS 

In the present paper, we demonstrated the existence of a new kind of spin 
current, connected with contraction and expansion of the Fermi surface. For 
one-dimensional current flow, it is larger than the traditional spin current based 
on translation of the Fermi surface. Each spin current is proportional to its 
own kind of spin imbalance. Finally, we introduced an equivalent dc electrical 
circuit describing the processes of electrical conduction, spin relaxation and 
spin- wave emission in a magnetic multilayer traversed by a dc electrical current 
perpendicular to layers. 

I thank Dr Ya. Bazaliy and Dr M. Tsoi for very useful discussions. 
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FIGURE CAPTIONS 

FIG. 1. a) Multilayer used for current-driven eperiments. Thick magnetic 
layer F\ with fixed magnetization prepares the spin of conduction electrons in a 
specific direction. On the other side of thin non-magnetic spacer N are the thin 
precessing magnetic layer Fi and the thick non-magnetic layer iVj. b) Definition 
of spin coordinates X2, in relation to the fixed atomic spin Si in F\ and 

to the precessing atomic spin S2 in F2. 

FIG. 2. a) The circles represent the spin-up and spin-down Fermi surfaces 
in N. A local electrical current density j x > causes translations by amounts 
(5fc+ < 0, 5k~ < of these Fermi surfaces along k x , with Afi < 0. b) Now, one 
Fermi surface is expanded and the other contracted. The radial distance 8k^ 
of the Fermi surfaces is proportional to the nonlocal spin imbalance Ajl < 0, 
according to Eq. (11). 

FIG. 3. Equivalent dc electrical circuit for conduction processes in the mul- 
tilayer of Fig. la. The upper horizontal resistors represent conduction in the 
spin-down band in the various layers, and the lower horizontal resistors in the 
spin-up band. The vertical resistors represent electron-spin relaxation between 
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the two bands. The dc electrical motor represents the interaction between elec- 
trons and spin waves in F 2 . 
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